A New Approach to Personal Income Distribution 

*Atushi Ishikawa a , ^Tadao Suzuki 6 and ^Masashi Tomoyose c 

a The Organization of Core Curriculum Studies, Kanazawa Gakuin University 

Kanazawa 920-1392, Japan 
6 Department of Management Information, Nanao Junior College 
Nanao 926-8570, Japan 
c Institute for Theoretical Physics, Kanazawa University 
Kanazawa 920-1192, Japan 

February 1, 2008 



Abstract 

The results of R 2 dynamical random surface model (2-dimensional quantum gravity with 
a R 2 term) are applied to explain the personal income distribution. A scale invariance exists 
if there is not the R 2 term in the action. The R 2 term provides a typical scale and breaks 
the scale invariance explicitly in the low and middle income range. A new distribution, 
Weibull distribution, is deduced from the action analytically in the low income range, and a 
consistent fitting is obtained in the whole income range. Also, we show that the lognormal 
distribution in the middle income range can be understood in this framework. 

1 Introduction 

The analysis of the personal income distribution is one of the important subjects in the econo- 
physics (p. The first step concerning this research is almost established. There are some 
empirical knowledge of it. It is believed that the distribution follows a power law in the high 
income range and follows a lognormal distribution in the low-middle income range. The purpose 
of this paper is to explain this behavior from a viewpoint of R 2 dynamical random surface (DRS) 
model. 

The distribution in the high income range is first discussed by V. Pareto in 1897 [gj. The 
probability density function of the personal income x is given by the following; 

p(x) = Cx- {1+a \ (1) 

where C is a normalization factor. Also a is the parameter which characterizes the profile of 
the income distribution, so-called Pareto index, and this power law behavior is known as Pareto 
law. 

On the other hand, it is indicated by R. Gibrat in 1931 that the probability density function 
in the low-middle range follows a lognormal distribution as follows Q; 



p ^ = 7?. — 9 exp 
xv lira* 
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where xq is a geometric mean value of x and a 2 is a geometric variance. 

In the early '80s, the above two knowledge are rediscovered 0. The income distributions in 
United States 1935-1936 are analyzed in these papers. Most of the data follows the lognormal 
distribution, while the higher 1% is governed by the power law. Recently, this fact is reconfirmed 
H by using the detailed Japanese data ||] . 

The personal income distribution has been an interesting object of study. What seems to 
be lacking is a model from which we can deduce the features in the whole income range. For 
example there are models employed a stochastic evolution equation. It is easy to derive a power 
law, however, it is hard to lead a lognormal distribution by using these models @. Because the 
basic idea of these models is that they contain no typical scale in the fundamental equations. We 
should comment an exception. In recent paper ||, the authors employed a concept of complex 
networks Q and a stochastic evolution equation with some interactions. As a consequence, 
they succeed to construct a power law and lognormal distribution in the whole range. There 
is, however, not an analytic consideration for creation of the lognormal distribution but only a 
phenomenological one by using computer simulations. 

Our viewpoint is simple as follows. There is no typical scale concerning to money in the high 
income range. The distribution of it, therefore, should obey a power law. On the contrary, there 
is a typical scale in the low-middle income range, and it causes the breaking of the power law 
distribution. As a toy model which can realize above two cases in the same framework, we employ 
R 2 DRS model which is one of the methods for studying 2-dimensional quantum gravity with a 
typical scale. 2-dimensional gravity without a R 2 term has no typical scale, unlike 4-dimensional 
gravity. Here R means the scalar curvature of a 2-dimensional surface. The dimension of R is 
(length) -2 , namely the R 2 term is dimensionfull. The R 2 term, therefore, provides a typical 
scale to 2-dimensional gravity. It is well known that this model has two ranges. A fractal and 
self-similar structure, i.e. a power law, is seen in the range where the typical scale does not affect 
the theory. On the other hand, the fractal breaks down in the other range where the scale has 
the sense. This peculiarity is similar to the profile of the personal income distribution. 

In order to confirm our view in the framework of R 2 DRS model, it is straight to adopt the 
money itself as the fundamental variable of theory not the money possessed by individuals. We 
have once forgotten the concept of an individual and pay attention only to the connection of 
the money in a certain unit. This connection forms a 2-dimensional surface in DRS model, and 
later the money of an individual is formed geometrically. 

The outline of this paper is as follows. Section |2| is a brief review of the formalism upon 
which our work is based. R 2 DRS is explained and MINBU analysis that is one of the methods 
for studying the property of random surface is described. The aspects of continuum quantum 
R 2 gravity in 2-dimensions and the two asymptotic forms of the partition function are presented 
in section ^. These two forms correspond to the distribution functions in the high and the 
low income ranges, respectively. The interpretation and relation between quantities in R 2 DRS 
model and ones in the personal income are discussed in section |j. In section ||, the practical 
fitting is shown. The personal income distributions of United States in the year 1935-1936 and 
Japan in the year 1997-2000 are employed for its purpose. Summary and discussion will be given 
in section |(| 



2 R 2 DRS Model and MINBU Analysis 

In this section, we review DRS and explain MINBU analysis which is one of the methods for 
studying DRS [10]. DRS model is known as dynamical triangulation (DT) |ll| in particle 



physics, and DRS with a R 2 term has a fractal and a non- fractal ranges ]l2). In DRS, we 
consider a 2-dimensional surface. This surface is constituted by attaching the sides of simplexes 
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(for instance equilateral triangles), and all possible constitutions are statistically averaged. 

The typical surface is shown in Fig. [j]. This surface has a self-similar structure. A simply 
connected region of area with one boundary is called a baby universe. In the dynamically 
triangulated surface, a minimum neck between baby universes is composed of three sides of 
triangles and can be identified uniquely. The distribution of minimum neck baby universe 
(MINBU) is, therefore, one of the important observable quantities in DRS (DT). In order to 
analyze it, we consider the case that a whole universe (area A) is divided into two baby universes 
(area A — B, B), paying attention to one of minimum necks (Fig. |2|). We denote that Z(A+ 1) 
is the fixed area A + 1 partition function, and that Z(A, I) is the partition function for surfaces 
(area A) with one boundary (length I). MINBU (for minimum boundary loop length / = 3) can 
be recognized as a surface from which one triangle is removed from a sphere. We can correlate 
these two partition functions, by counting ways which triangle is chosen and removed, as follows; 

Z(A,3)~(A + 1)Z(A + 1). (3) 

The partition functions of two MINBUs in Fig. || are given by 

Z(B,3) ~ (B + 1)Z(B + 1), (4) 
Z(A-B,3) ~ (A-B + 1)Z(A-B + 1). (5) 

We define, here, that ua(B) is the statistical average number of MINBUs of fixed area B on 
a closed 2-dimensional surface of area A. If the 2-dimensional surfaces once divided into two 
MINBUs are made to be one universe again, the following expression can be obtained; 

Z(B,3)Z(A-B,3) 
Z{A) • (6) 

This formula works well as long as we treat MINBUs. 

It is well known that DT is equivalent to the continuum theory [15|. If we take the continuum 
limit which brings every sides of triangles to the infinitesimal and the number of triangles to 
the infinity simultaneously keeping the area of a 2-dimensional surface constant, it is thought 
that this limit becomes equivalent to 2-dimensional gravity defined by the continuous variables. 
The explicit forms of the partition function Z and 714 (JB) is obtained by using the knowledge of 
2-dimensional continuum quantum gravity. In next section, we consider it. 

3 Continuum Quantum R 2 Gravity in 2-dimensions 

2-dimensional quantum gravity with a R 2 term is investigated in the continuum framework in 
ref . 1 13 ] . The partition function for fixed area A is given by 



where guuil^i v = 0, 1) is a metric, X l (i = 1, 2, ■ ■ ■ , d) is a scalar matter field and d is the number 
of scalar fields. Also R is the scalar curvature of the 2-dimensional surface. The 5-function fixes 
the area to A. The actions are defined as 

S(9) = ^/d^(^ 2 + 4, ), (8) 

S M (X % ;g) = i- J tfxy/gg^dpXidvXi, (9) 
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where (jlq and m are a bare cosmological constant and a coupling constant of length dimension 
— 1, respectively. The actions (||), (P) and the integration measure T>g,T>X in the partition 
function @ are invariant under 2-dimensional diffeomorphisms. The measure, therefore, should 
be divided by the volume of the diffeomorphisms, which is denoted by vol(Diff). The partition 
function ([?]) carried out by an appropriate gauge fixing has the exact expression. The partition 
function, however, contains infinite excited modes and is difficult to integrate the full functional 
completely. In Ref. O], instead of the full integration, the two asymptotic forms are shown: 



Z(A; m, jj) 



' 2tt 



const. A 7oo(d) -V 
const. A-wW^-V^ 



2tt 
m' 1 A 



for m 2 A — > oo, 
for m 2 A — * 0, 



where \i is a renormalized cosmological constant, and 

d - 25 - V(25 - d)(l - d) 



loo(d) 
7o (d, m 2 A) 



12 



+ 2, 



2 + 



(d- 12) 
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(10) 
(11) 



(12) 
(13) 



are string susceptibilities of 2-dimensional gravity. 

As for m 2 A — > oo, the asymptotic form of the partition function (|i~0| ) represents that the 
surface is fractal. The similar phenomenon can be seen in 2-dimensional gravity without the 
R 2 term [14], where no typical length scale exists. In this A ^> 1/m 2 range, even if the model 
contains the R 2 term, at an area scale much larger than 1/m 2 surfaces are smoothly fractal. 
On the other hand, as for m 2 A — > 0, the partition function (11) is highly suppressed by an 
exponential factor exp[ — ttt], hence, the fractal structure of 2-dimensional surface is broken. 
This phenomenon can be understood as follows. In the A <C 1/m 2 range, at an area scale much 
smaller than 1/m 2 surfaces are affected by the typical length scale, and they are not fractal. 



This is the model for which we were asking, however, the partition function (1C) cannot be 
made to correspond to Pareto law (jl|) directly. Because area A in Eq.(10) or (11) should be 
correspond to the whole money, on the other hand, x in Eq.(j^) or (|2|) is the individual income. 

In order to connect this model to the personal income, we should consider MINBU analysis 
in the previous section. By substituting the asymptotic forms (10), (11) to the partition function 
@, (H), We can find the following two asymptotic formulae; 



n A (B) 



const, 
const. 



17 B " 
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70-2 



exp 



2tt 



nv 



M) 



B 



for m 2 (l - f )£ 



oo, 



for m 2 (l - 4)B -> 0. 



(14) 
(15) 



From these expressions, we can recognize that the distribution of MINBUs obeys a power law 
in the limit m 2 (l — B/A)B — > oo, on the other hand, it does not follow any power law in the 
limit m 2 (l - B/A)B -> 0. 



4 The Map between R 2 DRS and the Personal Income 

In the previous section, we analyze how the 2-dimensional surface (the whole universe) is dis- 
tributed to MINBUs. Let us consider a relation between the quantities of R 2 DRS and those 
of the personal income. By comparison between Pareto law (|l|) and Eq. (14), we propose a 
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correspondence relation; 



(l-j)B & x, (16) 
n A (B) 



N 



& P(x), (17) 



where N is a normalization factor. We should note that MINBU is, in general, defined as 
the smaller universe when the whole universe is divided into two parts, therefore, B < A/2 
(x < A /A). From this correspondence, we can correlate the string susceptibility to Pareto index 
as follows; 

Too = 1 - a. (18) 

Since Pareto index a is reported as 1.5 in [l], ||, we should take the value of 7^0 as — |. 
This means that we should chose the corresponding 2-dimensional surface as a sphere without 
any matters (d = 0), and we can conclude that 70 = by using Eq. (|l3|). Two asymptotic 
distributions are, in the end, acquired; 

p(x) ~ const. x~2 for m?x — > 00, (19) 

27T 11 , 

for m 2 x -> 0. (20) 



~ const, x 2 exp 



m 2 x 



Here m is a typical area scale. We should notice that Eq. (p0|) is known as Weibull distribution. 

The asymptotic distribution (^TJ) is correspond to the low income one, and it is not a lognor- 
mal distribution which is offered by many researchers. This fact must be paid attention. The 



form of Eq. (20) obtained by our argument agrees with the practical data in the low income 



range. Particularly its fitting is better than the lognormal one in the low income range. 

5 The Data Fitting of the Income Distributions 

In this section, we attempt to fit the personal income distributions of Japan in several years 
1997-2000 |, §, and the United States in the year 1935-1936 Q. For this purpose, we define a 
cumulative probability as follows: 

P{x <)= dy p(y). (21) 

Jx 

Substituting Eqs. (|l9|) and (|20|) for this definition, we obtain below formulae: 



P(x <) ~ Cqo x 2 for m 2 x — » 00, (22) 



1 — Cq exp 



2vr 1 

m 2 x 



for m 2 x 0, (23) 



where Cqo and Co are normalization factors. 

Firstly, let us consider the distribution in the high income range. In our model, Pareto index 
a takes the value of 1.5. We find that this index is almost consistent with the result of |l], ||, ||] 
(see U.S. data in Fig. Q). Also we should notice that this value is little larger than that in Japan 
for the year 1997-2000 §] (see the Japanese data in Figs. | and |). 

Secondly, the profile in "the low" income range is considered. This paper gives the weight to 
analyze this distribution. It must be noted that "the low" does not mean the low and middle. 
This point of view is one of the main differences from previous works. The plots of the data and 
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the best fitting are also presented in Figs. || and ||. The parameter 27r/m 2 in the figures can be 
obtained as the slope in the graph where the horizontal axis is taken as 1/x and the vertical 
axis is taken as log(l — P(x <)) in the limit m 2 x — ► 0. 

In our model, the parameter m specifies the aspect of the breaking of a power law, and we 
are able to make a best fitting by changing the value of it. The best fit values of 27r/m 2 in Japan 
1997-2000 and in U.S. 1935-36 are given as 9.48, 10.8, 8.43, 7.93, and 2050, respectively (Figs, g 
and HJ) . We have employed the coefficient of determination r 2 for the parameter that indicates 
the appropriateness of practical fitting: 

_ m - (f)? 

-£(F,-<F»2- [M) 

Here F{ is the ith value of data, /j is the value of the fitting function and (F) is the mean value 
of F{. 

It remains an unsettled question what distribution function should be employed in the middle 
income range. We have to inquire, to some extent, into this subject. It is believed that the 
distribution function is a lognormal one in this range. It is difficult to deduce a lognormal 
distribution from the partition function (^) analytically. We employed, therefore, a computer 

simulation of DT with a R 2 term Pl^t ^ „^ 0] ^° investigate a profile in the middle 



3 t-l qi 

income range, i.e.x is not so high neither low. Here /3l is a coupling constant and qi is the 
number of links at the vertex in DRS. The result is shown in Fig. |5| as the data of the distribution 
of MINBUs which is not cumulated. The horizontal and the vertical axes are the logarithm 
of (1 — A/B)B and ua_{B), respectively. As the figure indicates, the profile in the middle 
range can be well approximated with a secondary curve, and this behavior deduces a lognormal 
distribution. It seems that we can obtain the whole personal income distribution from the 
partition function (Q). 



6 Summary and Discussion 

In this paper we have proposed the R 2 dynamical random surface (DRS) model as a toy model 
to explain the profile of the personal income distribution. The money in a certain unit have 
been identified with a certain unit area (simplex) of a 2-dimensional surface. The connection 
of the unit money have been corresponded to the attaching of simplexes' sides which forms a 
2-dimensional surface. On the basis of this correspondence, we have analyzed the distribution 
of the 2-dimensional surface statistically as the personal income distribution by using MINBU 
analysis. 

The point is that R 2 DRS model has a typical scale 1/m 2 , which is the coefficient of the 
R 2 term, and that the distribution of MINBUs exhibits two distinct phases. One is fractal in 
the large MINBU range where the typical scale does not affect the DRS, and the other is not 
fractal in the middle-small MINBU range where small MINBUs are suppressed by the effect of 
the R 2 term. On the other hand, this fact could be confirmed by the continuum theory which 
analytically presented the two asymptotic forms in these two region. 

Moreover, we have proposed that the personal income x should be corresponded to (1 — 
B/A)B, where A is the whole area of a 2-dimensional surface and B is the area of a MINBU. 
From this assumption, the power law distribution has re-obtained in the high income range 
similarly to Refs. [g, ||, [5], [7], [g. In addition, it has deduced directly from the partition function 
in our method. At the same time, we have identified that Pareto index a = 1 — 7oo 5 where 7oo is 
the parameter which characterizes the profile of the 2-dimensional surface in the limit x — ► oo, so 
called the string susceptibility. We have also concluded that we should chose the corresponding 
surface as a sphere without any matter fields (d = 0). 
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We have argued, in particular, the distribution functions in the low and the middle income 
range in detail. The new distribution function so-called Weibull function has been proposed 
in the low income limit. This function has well fitted to the data rather than the lognormal 
distribution. A method of least squares has been employed to evaluate a validity of the data 
fitting. Also we have employed the computer simulation to discussed the distribution function 
in the middle income range. Our whole results were displayed in Fig. ^ and Fig. |j. 

In our argument, we assume that the personal income x was identified with the area (1 — 
B/A)B. We try to explain this identification below. In MINBU analysis, a boundary of a 
MINBU is a minimum neck, where two universes are correlated through only three sides of 
simplexes each other (See Fig. ^). This means that the connection between universes is very 
weak at that point, therefore, it is natural to consider that the minimum neck separates the 
money which an individual possesses. 

We should notice that the area of a small MINBU is double counted as the area of bigger 
one by definition. This may be because why the income which an individual obtains should be 
added to the income of those who have obtained more income. The effect —B/A can be neglect 
when B is much smaller than A. We must, however, take into account of it when B is close to 
A/2. In our model, there is the huge mother universe from which all MINBUs are removed (See 
Fig. Q). It is hard to find the object corresponding to this mother universe in the real economics, 
however, the theoretical distribution obtained has well expressed the real data. There may be 
virtual heat bath of the money in the real world. These arguments are left behind as a future 
subject. 

In this paper, we have shown that the distribution of the personal income could be under- 
standable as the geometry of a 2-dimensional surface of a sphere without no matter field. This 
suggests that the technique developed for the analysis of spacetime structure can be applied also 
to the analysis of our social structure. In 2-dimensions, above geometry is simplest one and the 
geometry of a 2-dimensional surface with genus and matter fields is also well understood. There 
are many fractal objects in the econophysics, for instance the distribution of the company profit, 
and the extended geometry may suit these objects. We are pleased if this analysis becomes a 
first step for the research which considers our society geometrically. 
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Figure 1: Fractal 2-dim Random Surface 



Figure 2: Divided MINBs 



2000 Japan 




Weibull ( 27t/m 2 = 7.93 ) 
Lognormal ( a= 0.536, x = 3.81 ) 
Power law ( a = 1 .93 ) 



10° 10 1 10 2 

Income (millionyen ) 



Figure 3: Cumulative Probability for 2000 Japan 



9 



1997 Japan 







« 








• Income 

Weibull ( 2?tfm 2 = 9.48 ) 

Lognormal ( o= 0.499, x = 3.89 ) 

Power law ( a = 1 .96 ) 





10" 



10' 



Income (millionyen ) 
1999 Japan 



10 z 



• . . 

•v. 

».. 

• 




V 








• Income 
Weibull ( 2rfm 2 = 8.43 ) 

— Lognormal ( o= 0.547, x = 3.81 ) 

— Power law ( a = 1 .99 ) 





1998 Japan 



— io- 
-Q 
ctf 
-Q 
O 

> 

"■4—* 

CTJ 

£ 10~ 

O 



io- 



•* 

• 




< 


> 






• Income 

Weibull (2rfm 2 = 10.8) 

Lognormal ( o = 0.492, x = 3.97 ) 

Power law ( a = 1 .98 ) 





10" 



10' 



Income (millionyen ) 
1935-36 U.S. 



'3 

as 

_q , 

o 10 

Q_ 



~ io- 3 
as 

o 



10" 



10' 



10 3 10 4 10 5 

Income (dollars) 



Income (millionyen ) 

Figure 4: Cumulative Probability for 1997-1999 Japan and 1935-36 U.S. 





• 
• 








• 
























• Income 

Weibull { 2;r/m 2 = 2050 ) 

-- Lognormal ( o = 0.552, x = 1050 ) 
- Power law ( a = 1 .63 ) 





10" 



10 





log[(l-B/A)B] 

Figure 5: Lognormal fitting into the simulation of R 2 DRS model 



11 









Simulation ft = 100 
- Lognormal <r= 0.812 


















■Si, 

list 






















10° 10 1 10 2 10 3 



log[(l-B/A)B] 









Simulation ft = 100 
Lognormal <r= 0.812 


















■Si, 

list 






















10° 10 1 10 2 10 3 



log[(l-B/A)B] 



1 996 Japan 



V 

V 




< 








• Income 

Weibull(2rfm 2 = 9.14) 

Lognormal ( a = 0.512, x = 3.84 ) 

Power law ( a = 1 .79 ) 





10° 10 1 10 2 

Income (millionyen ) 



